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ABSTRACT 


Supply support of a Naval Air Rework Facility (NARF) 
should consider both on-site inventories of spare repair 
parts as well as back-up resupply from the local Naval Supply 
Center (NSC). This thesis presents a model for such a sys- 
tem for a limited time horizon. The decision variables are 
the number of units of an item to stock on-site and the length 
of time between deliveries once the on-site inventory is 
depleted. The determination of the optimal values of these 
variables required evaluation of the total expected variable 
costs for each given set of parameters. After identification 
of optimal values of both decision variables, a comparison 
between the minimum total expected costs of this model and 
an earlier model without on-site spares was conducted. The 
results suggest that the on-site spares model is preferable to 
one without spares. However, because the outcome of such a 
comparison is strongly dependent on the cost values assumed, 
additional analyses are needed before a general statement can 


be made. 
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I. INTRODUCTION 


A. BACKGROUND 

Mat Mwio, te nembcmartement Of Defense Material Distribution 
System Study recommended the consolidation of wholesale 
Supply support between the Naval Supply Centers (NSC) at 
Norfolk, San Diego, and Oakland, and their local Naval Air 
Stations. 

As part of a study of how to improve the local area 
material distribution at Oakland and San Diego, McMasters 
[Ref. l: pp. 1-6] addressed the problem of providing supply 
Support for the Naval Air Rework Facility in each air station. 
In that report, he suggested several support methods. One 
was to provide on-site inventories at the NARF. The advan- 
tages of such a plan would be quicker response time to the 
customer, reduced transportation costs, and reduced customer 
GeWawececeus. The costs of maintaining a separate inventory, 
however, would be a disadvantage. Another support method was 
through direct delivery by the NSC, with no on-site spares. A 
third possibility was a combination of the two methods; on- 
Site inventories at the NARF with direct deliveries from the 
NSC when demand exceeds on-site inventories. The optimum solu- 
Plone ommroolen Should be a trade off between customer 
needs, delivery costs, transportation costs, and costs which 


result from maintaining a separate inventory. 





In the process of modelling these alternatives, McMasters 
(Ref. 2: pp. 4-14] developed a model for determining the 
depth of repair parts to be stocked on-site using a time 
independent delay cost. This model also addressed costs of 
establishing an inventory on site, penalty costs for being 
out of stock at some time before the end of a fixed time 
period, and surplus costs for having units remaining at the 
end of the fixed time period. He also modelled three alterna- 
mive methods OF Dproviding direct deliveryesupport [Ref. l: 
pp. 7-42]. The costs included in these models were delivery 
costs, time dependent customer delay costs, and time indepen- 
dent customer delay costs. No attempt was made to combine 


the models. 


Bo PURPOSE 

This thesis will develop a model which combines both an 
inventory system for on-site spares and a direct delivery 
model. The direct delivery model to be used is the scheduled 
delivery model which provides for delivery of all demands at 
the end of every N periods if there is at least one demand 
during those N periods. The cost of a delay for this model 
will include a time dependent as well as time independent 
Sesusweeegoeresulting formulas will be subjected to parametric 
analysis for determination of optimum methods. The objective 
of the model will be to determine the on-site quantity which 


minimizes the total costs over N periods. 





Gee thsols ORGANIZATION 

Chapter I gives the background of the problem, the purpose, 
and the organization of this thesis. Chapter II establishes 
the probabilistic basis of the time dependent delay costs 
incurred by the NARF waiting for delivery from the NSC, given 
that there are N periods of time between deliveries. This 
result is conditioned by the probability that the on-site 
items are consumed before the end of N periods. Each element 
in the total cost equation is also examined, and the associated 
expected values are developed. Chapter III analyzes the model 
developed in Chapter II and seeks to identify the optimum 
number of on-site spares for various values of important 
parameters. Chapter IV compares the total cost results using 
no on-site spares and scheduled delivery as determined by 
Pavudcew| her. 3: pp. 18-35] to the optimal results of Chapter 
III. Chapter V summarizes the finding of Chapters III and IV 
and makes recommendations as the possible uses of the results 


and areas of further research. 





Pie SOpEL DEVELOPMENT 


This chapter develops the combination of an on-site 
inventory system and the scheduled direct delivery models of 


References 1 and 2. 


A. ASSUMPTIONS 
The work load schedule at a NARF is such that inductions 
of a component into overhaul occur at a specified constant 
rate during a quarter. As a consequence, the time between 
inductions is also a specified constant and can be used as a 
convenient measure of time. In this thesis it will be referred 
to as a time "period." 
Under scheduled deliveries without an on-site inventory 
system, the NSC's truck makes a delivery at the end of every 
N periods if there is at least one demand during the N periods. 
If, however, stock exists on-site then it can be used to fill 
demand and thus reduce the chance that the truck will need to 
make a delivery during N periods. 
In the development of the combined model the following 
assumptions will be made: 
1. The time between potential demands is the time between 
inductions of a component. 
2. At most one unit of a given repair part is required by 


each aircraft component undergoing rework at a NARF. 
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3. The probability that the component requires the repair 
pane 1S pa 

4. There are Y spare units of a given part on-site that 
will be consumed before requisitioning additional units 
from the supply center. 

5. Delivery will be made at the end of N periods if at 


least Y+l units are demanded during N periods. 


As an example of the process, consider Figure 2.1. There 
bess x INauctions per quarter, and Y = 1. The probability 
that a part will require replacement is 0.4. The first period 


induction requires a spare which is satisfied by the on-site 
Spare. The induction of periods 2 and 3 do not require spares. 
Another spare is required in period four and, since no stock is 
available, the part is requisitioned from the local NSC. An 
additional spare is required in period six creating two shortages 
GQuring the quarter. Since delivery is scheduled for the end 

of the quarter, delivery of both parts will be made to the NARF. 
The process will continue during the next quarter with the 


assumption of Y spares on-site. 


Q Z A. 6 


X= DEMAND) (O=NO_DEMAND! 


Figure 2.1 Example of Schedule Delivery with On-site 
Spareswith N = 6, Y = 1, and p = 0.1. 


dul 





Cost elements of the model will include Ct, the total 
round-trip costs of a delivery. These include the costs of 
a truck and a driver from the time it starts loading at the 
NSe@ £6 Ehe time it returns to the NSC from its delivery. 

There also will be special handling costs, Ch, and processing 
costs, Cp, incident to the expense of establishing a unit in 
inventory on-site at the NARF. 

Delay costs are costs incurred at the NARF as a conse- 
quence of not having a needed part the instant it is required. 
Two elements will compose the delay costs. One is S, which 
is the cost associated with putting a component aside. The 
second element is Cd, a time dependent delay cost, which will 
represent the delay cost per demanded unit per period. This 
element may include labor cost due to work stoppage, inventory 
holding costs, and cost associated with the non-availability of 


a repaired component to a fleet unit. 


Boe ee teRMiNISTIC DEMAND 
If demand from a customer occurs with certainty once 
every time period, the demand is deterministic. If a truck 
is dispatched every time a demand is received, the cost per 
unit delivered is Ct. If k demanded units are allowed to 
accumulate before delivery then the cost per unit delivered 
Pen eeamiine wt buck Capadeity iS nm units, then the delivery 
Ses OeL eames emimimized by walting until the truck is full. 
Pee mere wil) be eno delays as the on-site spares 


will satisfy all demand requirements of the system. However, 


iL 





if Y < N, and if we define k as (N-Y), then as k units accumu- 
Pate, the NARF experiences delay costs from lost production 
while awaiting parts. If we refer to the cost of delay of 

one time period as Cd, and if the NARF must wait until k 
demands are accumulated before the supply center makes a 


delivery, the total delay cost incurred by the NARF is 


Cait) 
2 

This formula is based on the deterministic demand of one 
unit per time period. If a unit is demanded during the first 
period after a delivery then there will be a delay of k-l 
periods until the arrival of the spare. Likewise, a demand 
Guede second time period after a delivery will have a delay 
See were pertods. SO the total delay time from delivery 
to delivery will be 


(k-1) + (k-2) + 2... + 1+ £0 KUke1) 


which can be rewritten, with the addition of the Cd tern, 


as 


C. RANDOM DEMAND 
Equation (1) assumes that a demand is made every period. 


litmomtismatetegeus EO a parti that 1s replaced in every component 


ibe 





undergoing an overhaul. if the repair part has a proba- 
Perey, Oe = 1, that it will be replaced in the component 
during overhaul, we then must consider a case where delay 

costs are a function of the probability of demand, as well 

as N and Y. Our goal is how to determine the expected costs 
per period associated with the system. 

The first stage of development will determine the expected 
total delay time given that there are Y repair items on-site. 
In order for a delay to occur, there must be Y+l units demanded 
before the end of the Nth time period. If we consider the 
process of replacement of a part of a component as an indepen- 
dent Bernoulli trial with probability p, then we must determine 
the distribution of the number of independent Bernoulli trials 
required for Y demands during N periods. This process is des- 
Seibed By Ehe negative binomial probability distribution. 

Let n be the number of trials (periods of demand), neces- 
Sary to observe Y demands. Clearly, the range of n is 
Pa wr) O. The negative binomial probability function 


ROG a Given Y 1S 


n-l, Y n-yY 
Perea —- (an): (1—p) (2) 
where 
me ey alee ee, 
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Bmore eseneatem, (2) represents the probability that it 
will take exactly n periods for a total demand of Y to occur. 
Seowecwamplc OLVAGuation (2), let us say there are a total 
Seeowcime periods, (N = G6), and the number of units on hand, 
Y, is one. The probability, p, of a demand during each time 
period is 0.4. The figure below illustrates the probability 
Pig@esele Ytne demand first occurs on each of the periods one 
through six. Since we have restricted the range of n to have 
an upper bound of N = 6, the sum of the probabilities will 


be less than one, 


0.1608 0.1920 0.1728 0.1382 0.1037 0.0744 


Brguae 22 ; 2 Distribution of Negative Binomial with 
mw — 17 N ="6,.-and p = 0.4. 


Given that it takes n periods for the Yth demand to 
occur, there will be exactly (N-n) time periods before a 
delivery is made if there is an additional demand before the 
end of the N periods. In determining the expected total 
delay cost for scheduled delivery with Y items on site where 
there are (N-n) periods remaining for a delay, it is useful 
to review the steps for obtaining the expected delay costs 


for N periods as presented in Reference l. 


ILS 





Bbelays are a function of the number of configurations 
that demands can take in N periods. The total number of 


configurations where exactly x demands occur is expressed 


by 


n — em (3) 


Pv-seOtal! number Gf configurations which can occur with at 


least one demand is 


To determine the expected total delay associated with the 
MeEcOmerourations given by (3), we first consider only those 
Somer ouraemons Naying exacely x demands where x > 1. The 


Beeopeabilatzy Of each such configuration is 
N- 
GZ) Sais (lea) ye (5) 


The number of configurations having a demand in period 


ie < }.-< News 
1 
1? : (6) 


heehee Cane theat em 15 independent of 3. Those 


eememeas ©Ceeurring in peziod 7 will have to wait until period 
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N for delivery and hence each must wait N-j periods. The 


total of all delays for those configurations having x demands 


is 


s =i 
ROS) eg) (N39) 
jel 
N 
N-1 
= (ay L (N-3) 
j=l 
TDix,w) = Sie (Noh) . (7) 


From (5) and (7) we can find the expected total delays over 


all x values: 


N 
Bren) v=.) TD(x,N)P(x,N) 
x=1 
N N(N-1) ,N-l Be N-xX 
=) x (Cj) oe ee) 
x=1 
Ni =i)ee N= 1 x N-x 
ETD(N) = 5 ) (Sia io SGlSs) (8) 
2 x~-l 
x=1 
By factoring the summed terms, 
N(N-1L) a N-l x-1 N-xX 
" ETD(N) rr rr (1-p) 
Z oe ill 
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The summation term is now equal to one and our expression 
reduces to 


erp(ny = NiX-lp | (9) 


Now, given that the last of the Y repair items is demanded 
in the nth time period, it easily follows that there are 
(N-n) time periods remaining in which demands may occur, the 
expected total delays can be expressed as 


ep(wen) = Geen) Liven) -2)p 


(10) 
However, since there is a probability associated with n time 
periods being required for the Yth demand to occur, Equation 
Peoemuceesbe Muliaplied by Ehat probability, given in (2), 
and the result summed over n to get the expected total delay 
over N time periods, given an on-site inventory of Y. The 


expected total delay cost is: 


~¥ , (N=n) ((N=n) -1)p 


5 jeans) 


fs _ . = es n 
ape Oe Cd) (je (1p) 
n=Y 


Upon examination of Equation (ll), we see that for a 
fixed N, the expected total delay cost is monotonically 
Gecreasing with an increasing Y. Furthermore, ETDC(N,Y) > 0 


as Y approaches (N-l). 
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Deeee oo lS INCURRED WITH THE ESTABLISHMENT OF Y ITEMS ON-SITE 

Reference 2 identified two costs, special processing and 
holding, associated with locating Y items on-site. Special 
peeeceocINg c@ses, Cp, are those paperwork and processing 
Geses Chargea when placing an item in an on-site store. 

If a quantity of Y items are placed into an on-site store the 
aeeal precessing Costs will be Cpy. 

The space required to store on-site units must be large 
pHougmecouaecemmodatc all Y units. In addition, the cost of 
that space can be expected to be constant even during periods 
when the number of items on-hand are less than the quantity Y. 


The total holding costs will therefore be Chy. 


Pee civ NDEPENDENT DELAY COSTS 

The cost per unit associated with putting a component 
aside when all Y spares are expended before delivery at time 
N will be denoted as S. This cost is assessed only at the 
time of the demand and is not therefore time dependent. It 
Can include the cost of placing a component in storage and 
documentation of the status of repair and requisitions. To 
find the associated expected delay cost we must first find 
the expected number of components that will suffer delays 
because of lack of repair parts when required. 

If x is the number of units demanded in N periods, then 
when x > ¥ the number of components set aside will be the 
Geeeee remem ye) The probability distribution that describes 


the total number of demands for a series of independent 


Ey, 





Bernoulli trials is the binomial distribution. Therefore the 


expected time independent delay costs are 


N-1 
S J) (x-¥) (2) p* (1-p)"™ . (12) 
xXx=Y+1 


Netce that the upper bound of the summation is N-1, for if a 
demand occurs in the Nth time period it is assumed to be 


filled immediately and no delay occurs. 


meee oURPLUS: COSTS 

Reference 2 also examines surplus costs. The unit cost 
for having a surplus of items in on-site inventory at the end 
of N time periods is kC where C is the unit cost and k is 
a factor which may be greater than 1.0. If x < Y, then the 
cost of surplus in N time periods is kC(Y-x). The expected 
total surplus cost may be found by the same method as the 


time independent delay cost. 
x vs 
own) () ee lpr as (13) 


Cele biveEryY COSTS 

When considering delivery cost, the cost associated with 
making a delivery at time N will be denoted by the term Ct. 
Since a delivery only can occur if the number of demands x 


exceeds the on-site spares Y in N time periods, the probability 
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gemxmexcececaing Y Will be the sum of the binomial probabilities 
for x from Y to N. Therefore the expected delivery cost will 


be: 


N 
oe ap) (14) 
x=Y+] ~* 


wee ee ne lp TOTAL COSTS WITH A FIXED N 
The expected total cost over a fixed total number of 
periods N is found by summing all of the expected cost 


elements described above. 


; 7 Spectaly dana Lang Surplus 
pew = Processing Costs’ ( Costs 
rn (Time Independent) (Delivery) 
Delay Costs Costs 
Time Dependent, | 
et Delay Costs re 
oe ae N-x 
Bier ee sep + eh)y + ke ) (1) p (1-p) 
x=0 
N 
N- 
+S ) ea) () p* (1-p)” * 
x=Y+1 
N 
et ) =) Dp leroy as = 
x=Y+1 
Sea n arene him) — 1) 
1 Gel Eee Gt) amen 
re 2 
gang 
(5G) 
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This chapter has developed a formula for expected total 
costs for a system which combines an inventory system for 
on-site spares and a direct delivery model. Equation (16) 
will next be analyzed in an attempt to determine the optimal 


values of the decision variables Y and N. 
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i ee CME ZALON: ANALYSIS 


Since N and Y can take on only discrete values, the use 
of finite differences is appropriate for determining their 
optimal values. However, since optimization formulas for N 
and Y based on finite differences were as complex as the 
Original cost equation (16), an APL program was written which 
numerically determined the expected total cost for a range of 
N and Y values. The program is included in Appendix A. The 
results were then plotted and optimal values were determined 
by examination of the graphical results. 

In order for this analysis to be comparable with Reference 
1 and Davdison's analysis in Reference 3, the following values 


were assumed for the cost terms: 


Time dependent delay cost (Cd) $50 per unit per period 
Time indepedendent delay cost (S) 920 per unit 
Seeceiaemandling cost (Ch) So Omecer Unit 

Specie emocessing cost (Cp) See OO oer unit 
SurpPuermcesc (kC) So 0m pem Unit 


Since the total number of time periods and probability of 
demand are likely to be fixed in practice, the analysis uses 
fixed values of N and p and varies Y. 

Figure 3,1) provides a look at the total cost of a scheduled 


delivery scheme with fixed N between from 10 and 50 and 
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10 20 30 40 50 
NUMBER OF ONSITE SPARES, Y 


Bigure ss. Motels Gost ascea punctlon Of On-site 
Spares with Varying Values of N. 
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p= 0.1. Although the total cost values are discrete, a curve 
is drawn through those points for clarity. The N values were 
chosen to illustrate the general shape of the curve for a 
range of N. This figure also illustrates the convexity of 
the curves when costs are presented on a linear scale. Subse- 
quent figures will use a log scale for total costs to facili- 
tate comparison between parameters and may not appear convex. 

Note that the number of on-site Spares, Y, was limited 
to a maximum value of N. This is consistent with the develop- 
ment of the model because if Y exceeded N then the number of 
on-site spares would always exceed total demand and this 
Seeiaeioneeleanrly would not lead to an Optimal solution. 
The cost curves illustrate this observation as the total 
cost approaches its maximum as the value of Y approaches N. 
As would be expected with a small probability of demand, 
the optimum value of Y is small relative to N. The time 
dependent delay costs are small with low demand and the 
surplus and special handling and processing costs will increase 
mech a large Y. 

Picmmeme-2 illustrates the components of the expected 
PertleecOscemeutye as described by Equation (16) for the case 
of N = 50. Figure 3.2 shows that the major components of the 
total cost curve are the time dependent delay cost, Cd, and 
the surplus cost, kC. For the given parameters, the delivery 
Sect Tele audmtme time imdependent delay cost, S, do not have 
a significant impact on total cost. As expected, both of these 


terms strictly decrease as Y increases. The special handling 
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and processing cost term, increases linearly with increasing 
Y but due to the value assigned to (Cp+cCh), the term has 
ieee Impact On total cost. 

Figure 3.3 presents total cost curves on a log scale for 
varying Y with a fixed N and three different values for the 
probability of demand, p. The first graph in 3.3 represents 
the same situation as Figure 3.1 but with costs presented 
with a log scale. As can be seen from the graphs, as the 
probability of demand increases so does the number of on- 
etewopares required for optimality for the same N. This is 
due to the increase in both delay cost terms when the demand 
exceeds Y, which is more likely with an increasing p. These 
costs will outweigh surplus costs for excess Y, which are 
less likely with increased probability of demand. 

Figure 3.4 illustrates how the optimal value of Y varies 
with N. The figure emphasizes the discrete values of N and 
Seetmumer oe Lae distinct break points for p = 0.1 become less 
pronounced as p increases, becoming nearly linear as p ap- 
Paodehecmie 0. Gwen p = 1.0, Y = N, both delay cost and 
surplus costs will be zero, and the only costs with a positive 
Pave eeeee the Ch and Gp terms, which are linear in Y. 
This result should apply whenever Ch and Cp are much less 
Bian Cdyss mend KC. 

Peeiremonomshiows ROwstme tocal Costs vary over a range Of 
probabilities Oo: demand for a fixed N and selected Y values. 


These graphs also show that for a given N value, systems 
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OPTIMUM NUMBER OF ON-SITE SPARES, Y 
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with low probability of demand require smaller on-site stores 
GO Minimize total cost, as was noted earlier. 

Figure 3.6 shows the relationship between optimal Y and 
the probability of demand for three values of N. Though p 
1s continuous, for the purposes of illustration, the proba- 
bilities of demand are varied from 0.1 to 1.0 with increments 
of 0.1. Again, as probability of demand approaches 1.0, 
the optimum value of Y approaches N. For the given set of 
cost parameters these plots show that optimal Y is approxi- 
mately equal to the expected demand pN. The relationship 
between Y and pN becomes more nearly linear with increasing 
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iV eco len oON) LTO A SCHEDULED DELIVERY 
So lenewETHOUrL ON-SITE SPARES 


Davidson in Reference 3 performed a parametric analysis 
of costs of a scheduled delivery system with no on-site 


Spares. The cost equation used by Davidson was developed by 


McMasters in Reference 1 and is of the form, 


ee [i atea ti Cd (N=1) p) --In{1-(1-p) 
n°  Toep 2 e 


(his) 


ECP(N) = ([ - ee ley) 
Equation (17) describes the expected costs per period and 
there are only two cost elements, delay cost and delivery 
cost. For a comparison between (16) and (17), an adjustment 


to (17) was made. Equation (17) was multiplied by N to 


obtain a total cost value over N and the time independent 


delay cost term from (16) was added. The result is Equation 
(dhe) 
N(N-=1) NON x N-x 
ETC(N) = cda[———F) +s JY ())p*(1-p) 
x=l1 * 
a N, _X N-xX 
meet 8) +(e) (i=p) : (18) 
x=l1 * 


The expected total costs of both (16) and (18) were 
numerically evaluated using the same values for Cd, S, and 


Ct as the examples presented in Chapter III. The Cp and 
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Ch terms in (16) were also the same as those used in Chapter 
IOC ES 

Table la presents the results of the comparison when the 
probability of demand is 0.1. It displays the total expected 
cost for the two models. Model 1 corresponds to Equation (16) 
with an optimal number of on-site spares, and Model 2 corres- 
ponds to Equation (18) for the same N values. Tables lb and 
He present the results for probability of demand of 0.5 and 
0.8, respectively. As Table 1 shows, the model employing an 
optimal number of on-site spares has a smaller expected total 
cost than a system not employing on-site spares and the 
savings provided by an on-site system can be significant. 
The total expected costs with no on-site spares start substan- 
tially higher for all three probability values, and increase 
faster with increasing N than they do for the on-site spares 


model. This difference can be explained by the impact of the 





Cd term. In the case of no on-site spares, the Cd term, 
CaN(N-1) p 
g. f 


increases at a geometric rate with an increasing N, whereas 
in the on-site model, cost savings are achieved by delaying 
the application of the Cd term. 

A key guestion in the future comparison of these two 
Models wilhiebe. “At what po2mt do the costs of the imple- 
mentation of an on-site system no longer make it preferable 


BOM SYSECM wlenOut On-Site Spares?” 
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ei, 1 


Gomearison of Model 1 and Model 2 with p = 0.1, 0.5, 0.8 


Table la 


Propaotmity © Demand = 6...) 


Number of Total Expected Total Expected 
Periods Costs Model 1 Costs Model 2 
10 MS Oe, VIL (Gia iss 
20 3052 65 1077.84 
30 432.70 235 Ono 
40 50 6rs5 4078.52 
50 Dod 6324.48 

Table lb 


Probability of Demand = 0.5 


Number of Total Expected Total Expected 
Periods Costs Model 1 Gosts Medel 2 
10 240.95 NR es ee 
20 37 150m! 5049.50 
30 485.14 IEA Sere 
40 59155 20 O00. 00 
50 693.14 Big? 5120 0 

Table lc 


Peemabilaeywor Demand = 0.8 


Number of TeEtal expected iigedl) Expected 
Pemroas Costs Model l Costs Model 2 
10 GS 2-4) one es 
20 Zoek S015 5.39 
30 S213) ial he O79 23 
40 SS) TS Ss BEI 39209 
50 463.7 49899 .48 
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Deotmeaminaeronuer Equation (16), we find that the 
spceeral processing ana handling costs, (Cp +Ch), and the time 
independent delay costs, S, can be examined to answer this 
emestaon. The condition of indifference between the two 


Mernmoas 1S described by equating (16) and (18), or, 


PRGererN) = “EPC(N) 


where Y* is the value of Y that minimizes total expected 
costs. 
BOuemeemen Gp + Ch) term and the S term will affect the 
value of Y* as they are varied. The effect of the (Cp +cCh) term 
is suggested as follows. Suppose that we represent Equation 


ro) by 


TVC = f, (Y) + (Cp+cCh)yY , (19) 
where f£(Y) represents all other cost elements of the equation 
Sxelucumam(cp -Ch). Lf Equation (19) were continuous in 
Y we could take the derivative with respect to Y and set the 


Besult eCaqual to zero, Or, 


af, (Y) 
oui 1 a 
a | epee ch = 0. (20) 





Rearranging terms gives, 
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Ces eg), 


1 = 
ee = (Cp + Ch) - C2) 
and Y* could be determined from Equation (21). Equation 


(21) clearly shows that a change in Cp +Ch will affect the 
wae Of Y*, 

ime determimation OF Ehe breakeven point for the (Cp +Ch) 
term can be found by fixing the value of all parameters with 
the exception of (Cp +cCh) and determining Y*. The expected 
meral cose Of (16) with Y* and N is compared to (18) with the 
same N. Tiesvabue of (Cp +Ch) 1s varied and Y is recomputed 
and again the total costs of the two equations are compared. 
The process continues until the value of (Cp +cCh) is found 
that makes the expected total cost of (16) and (18) equal. 

The effect of a change in S, the time independent delay 
Secroepewita respect tO Y* is Similar. Any change in S will 
change the shape of the total expected cost curve which in 
turn will affect Y*. The procedure described above can be 
meeG tO filma the breakeven point for S. 

Variation of the Ct and Cd values will have no effect on 
the preference of the model described by (16) over that des- 
pail cous memes) .) sUpOn examination of the Ct and the Cd terms 
Mie eucatommro), It 1S easy to see that the expected costs 
of these two component terms are maximized when Y = 0. Clearly 
at that point, the terms in Equation (16) are identical to 
their analog in (18). Therefore, a change in either Ct or 


Cd will not affect the preference of (16) over (18). 





V. CONCLUSION AND RECOMMENDATIONS 


ee | ~)6 CONCLUSION 

This thesis addressed the problem of supply support for a 
NARF by the local NSC. A model using a combination of 
scheduled deliveries from the NSC and on-site spares at the 
NARF was developed. The optimal results of the model were 
obtained for several values of the probability of demand. 
These results were also compared to an earlier model for 
scheduled delivery without on-site spares, developed in 
References 1 and 3. 

The most noteworthy point of this thesis is that the model 
that uses on-site spares was found to have significantly lower 
total expected costs than the model that does not use on-site 
spares. Although several of the cost values chosen were 
hypothetical since data for the delay costs does not exist, 
they do serve to provide a relative comparison. If the actual 
values of the cost parameters can be established, the analysis 
can be repeated easily to determine how much better performance 
can be obtained with a combined on-site spares/scheduled 
delivery model than with the pure scheduled delivery model, and 
the optimal number of on-site spares. 

Also of note is the utility of Figures 4 and 6 in the 
management of a supply support system. If such a set of 


figures were available for actual cost values and changes did 


SEs) 





owe wie to sprOod oi lity Of Gemand for an item, or the 
total number of time periods, the new optimal number of 


on-site spares could be readily determined from such figures. 


B. RECOMMENDATIONS 

McMasters in Reference 2 proposed two other delivery 
methods. One method assumed that a delivery is delayed until 
some fixed number of units of an item have been demanded. 
Delivery is then assumed to take place as soon as the last 
SemMam@emocclbsoume ame second method starts by counting time from 
when the first demand occurs after the truck has returned from 
the NARF and is ready for further deliveries. Delivery is 
made M-l periods after the first demand. These two models 
Should also be evaluated since they might provide lower costs 
than the current scheduled delivery model when combined with 
ene On-site stocking. 

However, comparison between the three combinations of 
on-site delivery when those inventories are depleted will 
require the evaluation of total variable costs per time period. 
Thus, a renewal argument will be needed. The basis for that 
argument has already been established by Reference 1. The 
first step of such an analysis is to relax the constraint 
iMoscoomimmeraoter El that a < N. dn fact, n has an infinite 
upper bound as Equation (2) has indicated. From Equation (2) 


the probability of no delivery in the first N periods is 


emcee) oo 
n=N 
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and the probability of no delay during the first delivery 


period is 
) : p(n;Y) . 22) 
N 


These wwo probability statements can be subdivided into 
a sequence of many periocs of length N. We can then consider 
=i PpOsoloi tienes Gf Y not being used up in N periods, 2N 
periods, etc., and extend the model of Chapter II to cover 
those mutually exclusive alternatives. 

Finally, an aspect which should be considered in these 
two Gelivery methods is when the on-site spares should be 
replenished as this event constitutes a renewal. Perhaps 
the next delivery after Y has been depleted should include Y 
units in addition to the demands which have occurred since 


eac Pigster mages Cepleted. 
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APPENDIX A 


APL PROGRAM FOR NUMERICALLY EVALUATING EXPECTED TOTAL COSTS 


7 GTEST 
[1] ARAARARAAAAAARARAAAAAARAAAARAARAARAAAAAARAAAAAARARAAAAAARARAAAAARARAAA 
[2] § THIS APL FUNCTICN COMPUTES THE EXPECTED TOTAL COST OF A 
(3) a COASIHED GN-SITE SPARES/SCHEDULED DELIVERY SYSTEM. THE TOTAL 
(4] & NUMBER OF FERIGDS N. IS INPUT BY TEE USER AND THE FUNCTION 
(5) a COMPUTES TOTAL COSTS FOR THE SYSTEM WITH THE NUABER OF ON-SITE 
Co] @ SPARES YARYIAG FRGH { TO 4. THE RESULTING ARSAY IS THEN EXAAIMED 
(7] 9 TO DETERMINE THE VALUE OF Y THAT YIELDS THE LEAST TOTAL COST. 
[8] # THE ARRSY OF Y VALUES WITH THEIR CORRESPONDING TOTAL COST ARE 
(9] #& FRIATED ALONG WITH THE VALUE OF Y THAT HIMIMIZES TOTAL COST. 
[19] AAAAAAAAAAAARARAAAAARAAAAAAAAARARAAARRARAARRARAAARARAAAAARAAARAAAAARAA 
Cif] a YARTAELE NAMES ASSIGNED TO PARAMETERS 


(7 oa heen alter PROEABILITY OF DEMAND 

() See Deeeccsa ss >= Tink GEPENDENT PELAY COST 

(ieee iaaa see "~> =" PROCESSING COST 

EGE tt 0) HOLDING COST 

BRT i aerate TIAE INDEFENDENT DELAY COST 

iGhmeiees,----—o=-=2 SURPLUS FEHALTY 

FO DELIVERY COST 

Civ] 

[25] AARAAAARAAAAAARAAAAARAAAARAHAARAAARAAAARAAAAARA AAA AAR AAAAAARAARAAAAAAA 
24] 


(22) aoe 
[23] ‘ENTE 
C24] Hed 
[25] 

[24] AARAARAAAATRAR AAAS AAR AR ARS AA RAR AAAAAAAAAAAAARAAAAA RAN RAA AAA AR AA AAA AAR 
(27] 

[23] a PARARKETER TNITIALIZSTION 
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